Dg-Meson as Quantitative Probe of Diffusion and Hadronization in Nuclear Collisions 
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The modifications of -Ds -meson spectra in ultrarelativistic heavy-ion collisions are identified as 
a quantitative probe of key properties of the hot nuclear medium. This is enabled by the unique 
valence-quark content of the Ds=cs which couples the well-known strangeness enhancement with 
the collective-flow pattern of primordially produced charm quarks. We employ a consistent strong- 
coupling treatment with hydrodynamic bulk evolution and nonperturbative T-matrix interactions for 
both heavy-quark diffusion and hadronization in the Quark-Gluon Plasma (QGP). A large enhance- 
ment of the Ds nuclear modification factor (Raa) at RHIC is predicted, with a remarkable maximum 
of ~1.5-1.8 at transverse momenta around 2 GeV/c. We show this to be a direct consequence of 
the strong coupling of the heavy quarks to the QGP and their hadronization via coalescence with 
strange quarks. We furthermore introduce the effects of diffusion in the hadronic phase and suggest 
that an increase of the _D-meson elliptic fiow compared to the Ds can disentangle the transport 
properties of hadronic and QGP liquids. 

PACS numbers: 25.75.-q 25.75.Dw 25.75.Nq 



> 

O 



The properties of nuclear matter at high tempera- 
tures and densities are under intense investigation. Nu- 
merical lattice simulations of Quantum Chromodynam- 
ics (QCD) predict the formation of a deconfined Quark- 
Gluon Plasma (QGP) at a transition temperature of 
Tc ~ 170 MeV P, 0. The QGP presumably fihed the 
early universe during the first few microseconds of its ex- 
istence. Ultra-relativistic heavy-ion collisions (URHICs) 
at BNL's Relativistic Heavy Ion Collider (RHIC) and at 
CERN's Large Hadron Collider (LHC) are aimed at pro- 
ducing the QGP and characterizing its properties. Thus 
far, the experimental results point to the creation of a 
color-opaque, strongly coupled and almost ideal liquid 
with initial temperatures well above Tc [1, However, 
this assessment remains qualitative at present 0| . A 
quantification of the physical properties of the strongly- 
coupled QGP in terms of its microscopic interactions calls 
for novel probes and observables. 

A special role in these efforts is played by heavy quarks 
(charm and bottom, Q=c,b). Their masses {rric^b — 
1.3, 4.5 GeV) are significantly larger than the typical tem- 
peratures of a few hundred MeV reached at RHIC and 
LHC. Therefore, their production is essentially restricted 
to primordial nucleon-nucleon collisions fi\, and their 
thermal relaxation during the subsequent evolution of the 
medium is delayed relative to light quarks by a factor of 
~ mq/T « 5-20 [sl-fl^. This renders the thermalization 
time of heavy quarks comparable to the typical lifetime 
of the QGP phase in URHICs, and thus makes them ex- 
cellent probes. Current heavy-quark (HQ) observables, 
mostly in terms of pertinent semi-leptonic electron-decay 
(e^) spectra [13l-[l5|, indicate a surprising degree of ther- 
malization and collectivity, especially through a large el- 
liptic flow [i^ ■ First data on _D-meson spectra con- 
firm this behavior (Tg} . 

The thermalization and collective flow imparted on 



heavy quarks by the ambient expanding QGP serves as 
a quantitative measure of their coupling to thermalized 
light quarks and gluons. Current experimental results 
call for an interaction strength that goes well beyond 
perturbative QCD (pQCD) with realistic strong coupling 
constants, as ~ 0.3 |17h24| . Calculations based on non- 
perturbative heavy-light quark interactions, followed by 
coalescence into D mesons [13, HO] , describe the observed 
modifications of spectra reasonably well, up to trans- 
verse momenta of ~ 5 GeV/c. However, since the 
spectra are a superposition of charm and bottom contri- 
butions, an enhanced discrimination power is desirable. 
In this Letter we argue that accurate measurements of 
Ds mesons considerably improve our ability to constrain 
the interactions of heavy fiavor in medium: by directly 
comparing D and Ds spectra and elliptic flow one can dis- 
entangle recombination effects in hadronization, as well 
as the charm coupling to the QGP and hadronic matter. 

One of our key points here is to test the dual role of 
the resonant HQ interactions with light quarks - figur- 
ing into both diffusion and hadronization - by utilizing 
the well-established enhancement of strangeness produc- 
tion in URHICs. The latter is amon g th e earliest sug- 
gested signatures of QGP formation [23. It manifests 
itself as a ~50% increase of strange-to-light hadron ra- 
tios in central nucleus-nucleus (AA) relative to pp col- 
lisions (e.g., in the K/t: ratio) [2y, [23|, and can be un- 
derstood as a strangeness equilibration within the statis- 
tical hadronization model [2^. Thus, if resonant HQ 
interactions and recombination with thermalized light 
quarks into D-mesons are key to explaining HQ observ- 
ables, heavy quarks must also couple to the equilibrated 
strangeness content of the QGP. In contrast to earlier 
predictions of inclusive D^ yields in URHICs [29, 30], 
we here perform a full dynamical treatment of charm 
diffusion and hadronization over a large range in trans- 



2 



verse momentum (pt), encompassing thermal and kinetic 
regimes. Most notably, we predict the nuclear modifica- 
tion factor (Raa) of the Ds to significantly exceed one, 
which would be the first of its kind for a meson at collider 
energies. The main, dynamical, information is, however, 
encoded in its pT-dependence, which allows us to scruti- 
nize both diffusion and hadronization effects, especially 
in comparison to D mesons. 

The Ds and D observables encode yet another aspect 
which has not been evaluated in the HQ context to date: 
the role of the hadronic phase. Following the common, 
empirically supported notion that multistrange hadrons 
(e.g., (j) and fl") decouple close to the hadronization tran- 
sition [m, [s^!, the same should apply to the Ds. Hence, 
after evaluating the effects of hadronic diffusion on the 
i?-meson, a comparison of D and Ds observables enables 
a quantitative assessment of the hadronic transport co- 
efficient. Since the HQ transport coefficient is believed 
to be closely related to the widely discussed viscosity- 
to-entropy ratio of QCD matter (rj/s) [l^, disentangling 
the former for hadronic and QGP phases would allow us 
to quantify the temperature dependence of ry/s. 

Let us start by describing our calculations of charm 
diffusion in medium. At the temperatures of interest 
in URHICs, T<400MeV, HQ kinetics can be approxi- 
mated by Brownian motion, with a thermal relaxation 
rate A(p, T) derived from a Fokker-Planck equation fsl- 
[13, [13) US Ell • The key input to calculate A are HQ scat- 
tering amplitudes; at low- and intermediate momenta, 
these are governed by elastic interactions of potential 
type. In the QGP, we employ thermodynamic T-matrix 
calculations [2Cl |33| with input potentials approximated 
by HQ internal energies computed in thermal lattice 
QCD (IQCD). The resulting T-matrices exhibit resonant 
states close to threshold in mesonic and diquark channels 
close to Tc- They accelerate thermal HQ relaxation by 
a factor of ^3-5 over leading-order pQCD (s^]. In the 
hadronic phase, the thermal relaxation rate of D mesons 
is obtained 34] from elastic scattering amplitudes based 
on effective hadronic Lagrangians constrained by vacuum 
spectroscopy. Figure [T] summarizes our input for the 
spatial diffusion coefficient of charm in matter, defined 
by the Einstein relation as Vs = T/[mc,DA{p = 0,T)]. 
Coming from high T we observe a steady decrease in the 
QGP (due to an increased interaction strength in the 
IQCD potentials), followed by an increase in hadronic 
matter. The calculations support a continuous evolu- 
tion through the transition region, with a minimum of 
Vs ~ (4— 5)/27rr around Tc (translating into an esti- 
mate oi rj/s ^ (2— 4)/47r). Also note that our results 
in the QGP are comparable to recent quenched IQCD 
estimates [H, H^, but far below the pQCD values. 

The Brownian motion of charm is implemented into 
URHICs via relativistic Langevin simulations, with the 
space-time evolution of the medium approximated by 
boost-invariant ideal h ydr odynamics. To this end we em- 
ploy our recent tune [S^ of the AZHYDRO code [38| . 
optimized to describe bulk and multistrange hadron 
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FIG. 1: (Color online) The spatial diffusion coefficients "Ds (in 
units of the thermal wavelength, l/(27rr)) vs. temperature (in 
units of Tc) for charm quarks using T-matrix interactions in 
the QGP (red band) and D mesons using effective lagrangians 
in hadronic matter (blue solid line), compared to pQCD (pur- 
ple dashed line) and quenched IQCD (data points) [ssl, [3q|. 



spectra and elliptic flow in Au-Au collisions RHIC. It 
utilizes a state-of-the-art IQCD equation of state 0, 
[39} with pseudo-critical deconfinement temperature of 
Tc=170MeV, and a subsequent hadron-resonance-gas 
phase with chemical freezeout at Tch=160MeV to ac- 
count for the observed hadron ratios. A compact ini- 
tial spatial profile with pre-equilibrium flow permits a 
simultaneous fit of multistrange- and bulk-hadron data 
at chemical and thermal (Tfo~110MeV) freezeout, re- 
spectively. The initial HQ distributions are taken from a 
PYTHIA tune to spectra in pp and dAu collisions fl^l ■ 
The Cronin effect in nuclear collisions is accounted for 
via a Gaussian transverse-momentum broadening with 
{k\) = 0.6 GeV^ estimated from recent PHENIX e± 
spectra in dAu [iO]. 

After diffusion through the QGP charm-quark distri- 
butions are converted into charmed hadrons. We ac- 
complish this by applying resonance recombination (4]| 
with thermal light and strange quarks into D and Dg 
mesons on the hydro hypersurface at Remain- 
ing c-quarks are treated with i5-function fragmentation 
(as used in the fits to pp data). For a reliable coales- 
cence dynamics at low and intermediate pt, it is cru- 
cial that the formulation of the resonance recombina- 
tion model (RRM) via a Boltzmann equation [4l| yields 
the long-time limit of thermal equilibrium. We have 
verified this in the present case with the full space- 
momentum correlations as given by the hydrodynamic 
flow field jl^. The coalescence probabilities are esti- 
mated via Pcoai(p) — ATrosrc°''(p), with the charm-quark 
reaction rate, V^^ij)) (as given by the heavy- light T ma- 
trix), and a time duration Axres characterizing one gen- 



eration of D and Ds resonance formation [42|. With 
r;;'="(0)~0.2GeV [43I and Arros«lfm/c, we assume a re- 
combination probability of one at vanishing charm-quark 
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momentum, decreasing thereafter as determined by the 
dynamics of the RRM expression [41]. The latter is evalu- 
ated with mc=1.7GeV, mu.d=Q-3GeV, ms=0.4GeV and 
mr,=2.1GeV, m^?, =2.2 GeV with Tdm, =0.2GeV, ap- 
proximately representing the in-medium values of the T- 
matrices [33| during the hadronization window. Since 
HQ resonant scattering is underlying both diffusion and 
hadronization interactions, there is, in principle, some 
overlap between the two (this does not apply to the 
non-resonant parts of the interactions). To character- 
ize this uncertainty, we will study a scenario where dif- 
fusion interactions in the QGP are completely switched 
off for about 1 fm/c prior to Tc, corresponding to a tem- 
perature window of 180-170 MeV. The Langevin simu- 
lation resumes with hadronic diffusion of the combined 
coalescence-f fragmentation distribution for D-mesons for 
T<Tc until hydrodynamic freezeout at Tto=110MeV, 
while the ZJ^-meson distribution is frozen at Tc- 

It remains to determine absolute magnitude of the co- 
alescence contribution to the D and Dg yields in AA 
collisions. In pp collisions we assume fragmentation only 
with hadronization fractions from recent PYTHIA sim- 
ulations [H], i.e., D/c=82% and Ds/c^n%, includ- 
ing feed-down from excited states (here, L)=L)+-t-Z?°). 
This gives Dg/ D— 0.134 in pp, in line with CDF data 
m pp(yi=1.96TeV) MM\ 

and the value used in a re- 
cent PHENIX analysis [l5[. Since our coalescence contri- 
bution is evaluated using thermalized light and strange 
quarks within RRM, the logical choice for the pertinent 
Ds/D ratio are thermal weights, which we adopt from 
the statistical hadronization model (including feeddown 
from excited states}J[47|, with an additional strangeness 
fugacity, 7s =0.85 [291, for consistency with the hadronic 
equation of state in our hydro evolution (37| . Then, 
upon combining coalescence and fragmentation with the 
probabilities elaborated above, we obtain D/c—75% and 
Ds/c=15%, or Ds/D=0.20, for Au-|-Au collisions at 
^snn=200 GeV, hence obtaining an enhancement of the 
Ds/D ratio of ~50% over the value in pp. 

Figure [2] summarizes our D and Dg meson spectra in 
semi-central Au-|-Au at RHIC relative to pp collisions 
in terms of the nuclear modification factor, RAAiPr) = 
{dN^^/dpT)/{NcoiidNPP/dpT) (iVcoii: number of binary 
NN collisions in AA), and elliptic-flow coefficient, V2{pt) 
(the second harmonic of the azimuthal-angle depen- 
dence). Both the D and Dg Raa (upper panel) exhibit 
a maximum around pr— 2-3 GeV, induced by the trans- 
verse flow picked up from the expanding medium. Cur- 
rent STAR data are consistent with our Z?-meson result, 
but we predict the maximum to be more pronounced for 
the Dg (reaching beyond 1.5) due to c-quark coalescence 
with the enhanced strangeness in Au-|-Au. To further il- 
lustrate this effect, we also plot the result for c-quarks at 
the end of the QGP phase, which would directly represent 
D- and Dj-spectra if coalescence were absent and only 
(5- function fragmentation applied (as in pp). One clearly 
recognizes the important effect of coalescence, which only 
ceases above px—Si GeV, where fragmentation takes over 
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FIG. 2: (Color online) Our results for the nuclear modifica- 
tion factor (upper panel) and elliptic flow (lower panel) of Ds 
(red bands) and D mesons (green dash-dotted lines) in semi- 
central Au+Au collisions at RHIC. We also show the result 
for charm quarks at Tc (purple dashed lines), the equilibrium 
limit for Ds mesons in the hydrodynamic medium at Tc (blue 
dash-double-dotted line) and preliminary STAR data for 
the D-meson Raa in 0-80% Au+Au. In the upper panel, the 
red uncertainty band is governed by inclusion or omission of 
a Cronin effect in the initial charm spectra, while in the lower 
panel it is controlled by including or neglecting diffusion ef- 
fects in the hadronization window. 



and the D, Dg and c-quark Raa merge. While the c- 
quark spectra are not observable, the D and Dg ones 
are, so that their difference gives a quantitative measure 
of the coalescence effect. It turns out that hadronic diffu- 
sion does not significantly affect the Z3-meson Raa (due 
to a compensation of a decreasing temperature and an 
increasing flow of the medium). 

The elliptic flow of particle spectra is known to be an 
excellent measure of the medium's collectivity due to hy- 
drodynamic flow in non-central AA collisions (induced by 
the "almond-shaped" initial nuclear overlap zone). In our 
calculations, the diffusion in the QGP imparts an appre- 
ciable V2 on the charm quarks of up to ^4.5%, cf. lower 
panel in Fig. [21 Coalescence with thermal quarks ampli- 
fies this value by up to 50%, for both D and Dg mesons. 
However, while the Dg spectra freeze out after hadroniza- 



4 



tion, the D coupling to the hadronic medium, which in- 
herits the fuh eUiptic flow from the QGP expansion [37| . 
further augments V2 by up to 30%. We therefore sug- 
gest the U2-sphtting between D and Ds, in the spirit of 
the early multistrange freezeout in the underlying hydro 
evolution, as a promising measure of the transport prop- 
erties of the hadronic phase. 

In summary, we have argued that measurements of D 
vs. ZJg-meson Raa and V2 in URHICs provide power- 
ful constraints on heavy-flavor diffusion and hadroniza- 
tion. We have made predictions for these observables 
employing a self-consistent framework where the con- 
cept of a strongly coupled QGP is implemented in both 
macro- and microscopic components of the calculation: 
a hydrodynamic medium evolution, quantitatively tuned 
to bulk- and multistrange-hadron observables, has been 
combined with nonperturbative charm-diffusion coeffi- 
cients in the QGP which are compatible with currently 



available IQCD results. The diffusion of D mesons in the 
hadronic phase has been implemented for the first time, 
while Ds mesons are frozen out at Tc. A remarkable en- 
hancement of the Ds-iiieson Raa well above one emerges 
as a result of a strong charm-quark coupling to the QGP 
and subsequent recombination with equilibrated strange 
quarks. The latter can be directly tested by comparing 
Ds- and ZJ-meson Raa- The D meson picks up signif- 
icant additional V2 from the hadronic phase, which can 
be quantified by comparing to D<,-meson V2 for which 
hadronic diffusion effects are absent. This picture should 
persist at LHC and can be directly carried over to the 
bottom sector using B and Bs mesons. 

Acknowledgments. — This work was supported by the 
U.S. National Science Foundation (NSF) through CA- 
REER grant PHY-0847538 and grant PHY-0969394, by 
the A. -V. -Humboldt Foundation, and by the JET Collab- 
oration and DOE grant DE-FG02-10ER41682. 



Y. Aoki et al, Nature 443, 675 (2006). [23 

A. Bazavov et al, Phys. Rev. D85, 054503 (2012). [24 
J. Adams et al. [STAR Collaboration], Nucl. Phys. A757, 

102 (2005); K. Adcox et al. [PHENIX Collaboration], [25 
Nucl. Phys. A757, 184 (2005). 
K. Aamodt et al. [ALICE Collaboration], Phys. Rev. [26 
Lett. 105, 252302 (2010); G. Aad et al. [ATLAS Col- 
laboration], Phys. Rev. Lett. 105, 252303 (2010). [27 
M. Gyulassy and L. McLerran, Nucl. Phys. A750, 30 
(2005). [28 
E. Shuryak, Prog. Part. Nucl. Phys. 62, 48 (2009). 
S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 94, 082301 (2005). [29 

B. Svetitsky, Phys. Rev. D37, 2484 (1988). 
H. van Hees and R. Rapp, Phys. Rev. C71, 034907 [30 
(2005). 

CD. Moore and D. Teaney, Phys. Rev. C71, 064904 [31 
(2005). 

M.G. Mustafa, Phys. Rev. C72, 014905 (2005). [32 
R. Rapp and H. van Hees, in R.C. Hwa, X.-N. Wang 
(Eds.), Quark Gluon Plasma 4, World Scientific, 111 [33 
arXiv:0903.1096 [hep-ph]]. [34 
B.I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 
98, 192301 (2007); [Erratum-ibid. 106, 159902 (2011)]. [35 
A. Adare et al. [PHENIX Collaboration], Phys. Rev. [36 
Lett. 98, 172301 (2007). 

A. Adare et al. [PHENIX Collaboration], Phys. Rev. [37 
C84, 044905 (2011). 

Y. F. Zhang, J. Phys. G38, 124142 (2011). [38 
H. van Hees, V. Greco and R. Rapp, Phys. Rev. C73, 
034913 (2006). 

S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, 
Nucl. Phys. A784, 426 (2007). [39 
A. Adil and I. Vitev, Phys. Lett. B649, 139 (2007). [40 
H. van Hees, M. Mannarelli, V. Greco and R. Rapp, Phys. [41 
Rev. Lett. 100, 192301 (2008). [42 
P.B. Gossiaux and J. Aichelin, Phys. Rev. C78, 014904 [43 
(2008). [44 
S.K. Das, J.-e Alam and P. Mohanty, Phys. Rev. C80, 
054916 (2009). [45] 



W.M. Alberico et al, Eur. Phys. J. C71, 1666 (2011). 
J. Uphoff, O. Fochler, Z. Xu and C. Greiner, Phys. Rev. 
C84, 024908 (2011). 

P. Koch, B. MuUer and J. Rafelski, Phys. Rept. 142, 167 
(1986). 

B.I. Abelev et al. [STAR Collaboration], Phys. Rev. C79, 
034909 (2009). 

S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. 
C69, 034909 (2004). 

P. Braun-Munzinger, K. Redlich and J. Stachel, in 
R.C. Hwa et al. (eds.), Quark gluon plasma 491 [arXiv: 
|nucl-th/0304013| . 

A. Andronic, P. Braun-Munzinger, K. Redlich and 
J. Stachel, Phys. Lett. B571, 36 (2003). 

I. Kuznetsova and J. Rafelski, Eur. Phys. J. C51, 113 
(2007). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C79, 
064903 (2009). 

M. He, R.J. Fries and R. Rapp, Phys. Rev. C82, 034907 

(2010) . 

F. Riek and R. Rapp, Phys. Rev. C82, 035201 (2010). 
M. He, R.J. Fries and R. Rapp, Phys. Lett. B701, 445 

(2011) . 

H.T. Ding et al. J. Phys. G38, 124070 (2011). 

D. Banerjee, S. Datta, R. Gavai and P. Majumdar, Phys. 

Rev. D85, 014510 (2012). 

M. He, R.J. Fries and R. Rapp, Phys. Rev. C85, 044911 

(2012) . 

P.F. Kolb and U. W. Heinz, in R. C. Hwa 

et al. (eds.), Quark Gluon Plasma, 634 

[arXiv:nucl- th/0305084 ;P.F. Kolb and R. Rapp, 

Phys. Rev. C67, 044903 (2003). 

S. Borsanyi et al, JHEP 1011, 077 (2010). 

J. Durham, private communication (2012). 

L. Ravagh and R. Rapp, Phys . Lett . B655, 126 (2007). 

M. He, R.J. Fries and R. Rapp,|arX^:1106.6006 [nucl-th]. 

F. Riek and R. Rapp, New J. Phys. 13, 045007 (2011). 

G. Martinez-Garcia, S. Gadrat and P. Crochet, Phys. 
Lett. B663, 55 (2008). 

D. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 



5 



91, 241804 (2003). J. Stachel, Phys. Lett. B659, 149 (2008). 

[46] Y. Akiba, private communication (2012). 
[47] A. Andronic, P. Braun-Munzinger, K. Redlich and 



